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We develop a simple continuum model to analyze the vibrational modes of a nanomechanical 
multi-element structure. In this model, arrays of sub-micron cantilevers located symmetrically on 
both sides of the central clamped-clamped nanobeam are replaced by a continuum. In this approach, 
the punctual shear forces exerted by the cantilevers on the central beam are smoothed out and the 
equations of motion of the structure become exactly solvable. Our analytical results capture the main 
features of the vibrational modes observed both numerically and experimentally. Furthermore, using 
a perturbative approach to treat the nonlinear dynamics of the structure, we establish its frequency- 
amplitude response and analyze the mechanism of anharmonic coupling between two specific widely 
spaced modes of the resonator. 
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I. INTRODUCTION 

Nanomechanical resonators have been used to inves- 
tigate fundamental physics problems in a wide range of 
research areas. These include quantum measurement and 
quantum computation [TJ |5J |3J HJ [5] , ultra-sensitive force 
and mass detection [6], single spin detection [7], grav- 
itational wave detection [H [5] and other fundamental 
phenomena [5] |5J EH E]. The central reason for the 
increasing interest and activity in nanomechanical sys- 
tems for quantum studies is straightforward: fast dynam- 
ics of nanomechanical systems enable the investigation of 
the yet-unexplored corner of the parameter space where 
ultra-high frequency resonance modes cooled to near-zero 
temperatures are expected to manifest quantum mechan- 
ical corrections to the classical motion [TU [T3J Q3] . 

Likewise, from a technological standpoint, compact 
size, robustness, and high frequencies of nanomechani- 
cal resonators have resulted in numerous proposals for 
applications where technologies involving nanomechani- 
cal systems can offer alternative solutions to the existing 
electronic circuits and off-chip devices: frequency selec- 
tive oscillators [15], passive filters [16] . memory elements 
[T7] and spintronics devices [15] , 

Despite the broad areas of interest in nanomechani- 
cal systems, comprehensive studies of elasticity and me- 
chanics of nanomechanical structures are yet to be done 
beyond the simple single-element structures such as can- 
tilevers and clamped beams. The main difficulty has been 
the extension of elasticity theory to multi-clement struc- 
tures for closed- form analytical solutions. As described 
in the following, most of the eventual applications of 
nanomechanical systems will include complex structure 
design. Therefore, analytical understanding of resonant 
modes and other dynamical aspects such as nonlinear 
mode coupling in complex multi-clement structures will 
be of tremendous importance. Here, we briefly describe a 
few areas of interest in which nanomechanical structures 
in the gigahertz range will play a key role. 



A. Entry into the Quantum Regime 

Experimental access to the quantum realm is crudely 
defined as the regime in which the quantum of energy hf 
in a resonant mode with frequency / is larger than the 
thermal energy ksT of the environment. The motiva- 
tion behind this crude definition of the quantum regime 
is simple. The motion of a damped mechanical system 
can be described by a harmonic oscillator potential. In 
the quantum regime, the harmonic oscillator potential 
energy levels are discrete. In order to observe effects of 
discrete energy levels, smearing by thermal energy — due 
to finite temperature of the mechanical system — must be 
small compared to the energy level spacing, hf. The con- 
dition hf > kgT gives physically relevant parameters: 
a nanomechanical structure with a normal mode reso- 
nance frequency at 1 GHz will enter the quantum regime 
below a temperature T = (/i/fee)/ = 48 mK. Since 
typical dilution cryostats have a base temperature of 10 
mK, nanomechanical structures with frequencies above 
1 GHz can enable experimental access to the quantum 
regime. The experimental challenge is then to fabricate 
structures capable of high gigahertz-range resonance fre- 
quencies, and to measure them at low millikelvin-range 
temperatures. Since the resonance frequency increases 
with decreasing size of the system, one or many of the 
critical dimensions of the gigahertz-frequency oscillators 
must be in the sub-micron or nano scale. 



B. Mesoscopic Elasticity: Multiscale Modeling 

For large mechanical structures, continuum mechanics 
of elasticity theory provides the appropriate system re- 
sponse, both static and dynamic. As system size shrinks 
down to the submicron or nano scale, the elastic behavior 
starts to become atomistic rather than continuous, and it 
gives rise to a host of anomalous behaviors [19J. These in- 
clude enhancement of influence of surface defects, novel 
dissipation mechanisms, reduction of the stiffness con- 
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stant, and statistical fluctuation effects. Nanomechan- 
ical resonators therefore need to be modeled by atom- 
istic simulation in order to capture the essential aspects 
of their elastic properties arising due to mesoscopic size. 
However, atomistic simulation of these structures con- 
taining about 100 million atoms or more becomes com- 
putationally intensive. Further problems arise due to the 
existence of multiple length scales. Modeling of struc- 
tures too small for continuum mechanics and too large 
for atomistic molecular dynamics simulation requires a 
proper understanding of the coupling of length scales. 

For a comprehensive study of nanomechanical sys- 
tems, the structures need to be properly characterized 
for their elastic behavior. Specifically, correlation of sim- 
ulated modeshapes and experimentally-measured modes 
will require the knowledge of relevant static and dynamic 
parameters. Therefore, analytical studies of complex 
nanomechanical systems will not only provide a strong 
motivation for new approaches to model materials beyond 
the current limit of computing capacity, it will augment 
the ongoing work on the multiscale modeling of fracture 
dynamics and nanomaterials. 



C. Nanomechanical Device Applications in the 
Gigahertz Range 

The driving force behind MEMS (Micro-Electro- 
Mechanical Systems), of which NEMS (Nano-Electro- 
Mechanical Systems) or nanomechanics is the natural 
extension, has been the portfolio of MEMS applications 
in optical communication (routers, switches, repeaters), 
passive devices in cell-phone industry (filters, accelerom- 
eters, capacitors and inductors for integrated chip de- 
sign), and sensor technologies in chemical, biomedical, 
and electrical solutions. Nanomechanical structures are 
faster (because of higher resonance frequencies) , therefore 
the natural area of growth for these applications includes 
technologies for multi- function nanomechanical devices: 
ultra-sensitive force sensors for the detection of funda- 
mental forces and biomolecular forces. The technolo- 
gies needed for observing quantum effects will have obvi- 
ous use in these and other gigahertz-range applications. 
In fact, the range of frequencies discussed in this paper 
matches with the frequency bands for communication in 
a number of consumer devices: Cellular (0.4 GHz, 0.85 
GHz, 1.9 GHz), WiFi and Bluetooth (2.4 GHz), Satellite 
radio (~ 2.3 GHz), and GPS receivers (L1-L3: 1.2-1.6 
GHz). 

In basic research, hybrid nanomechanical devices in 
the gigahertz range can prove to be crucial in a num- 
ber emerging fields. For example, there has been a recent 
proposal for spintronics based entirely on nanomechan- 
ical torque associated with the electron spin [TS]. This 
device is an example of multi-function nanomechanical 
device, which, in the gigahertz range, can provide access 
to control and manipulation of electron spin at the spin 
coherence and spin relaxation time scales of nanoseconds. 

Since the structure size is reduced for increasing res- 



onant mode frequencies — and hence the operation speed 
of the device, it is possible that within a few years certain 
applications with nanomechanical devices will emerge 
which will have frequencies high enough so that at the 
relevant operating temperatures the nanomechanical ele- 
ment will be quantum mechanical. A specific application 
of gigahertz-range nanomechanical oscillators involves 
space communication devices (0.5-18 GHz). As passive 
devices, nanomechanical structures in certain space com- 
munication applications are expected to remain at low 
temperatures, perhaps near the quantum regime. This 
leads to the obvious realization that further shrinkage 
of micro- or nanomechanical systems with corresponding 
high gigahertz-range frequencies will inevitably require 
new paradigms. 



D. Approaches to Gigahertz Resonance 
Frequencies 

Central to many of the aforementioned applications of 
NEMS resonators is the attainment high natural frequen- 
cies of motion up to and beyond the 1-GHz mark. Among 
the various solutions proposed for achieving ultra-high 
frequencies in nanomechanical resonators such as the use 
of high stiffness materials [221 EI] and bulk mode geome- 
tries |22j , our approach of coupling mechanical elements 
to enhance high order resonant modes of the resonator 
structure has been shown to offer a number of advantages 
for performance and detection ease [23 . The coupling 
of additional degrees of freedom to a 10-/im-long simple 
nanobeam structure can produce significant modification 
of the high order resonance spectrum, with enhanced am- 
plitude and quality factors of selected vibration modes 
extending well beyond the 1 GHz frequency range. 



E. Organization of the Paper 

Here, we present two analytical models that yield 
closed-form solutions describing the dynamics of the 
coupled-beam resonator, dubbed the antenna structure. 
This structure is a prototype of a class of two-element 
structures that can be envisioned in specific applications. 
The continuum model allows for a clear comparison of 
the modal shapes and spectrum with full finite clement 
analysis of the structure. We further investigate nonlin- 
ear modal coupling between widely spaced modes of the 
structure using perturbation theory techniques. 

In Section II, we describe a discrete model of the 
antenna-like multi-element structure, and we obtain ex- 
act solution for the vibrational modes. In Section III, we 
extend this analysis to a continuum approximation where 
the periodic and punctual force densities are smoothed 
out . In this continuum model, we describe the frequency 
spectrum and the corresponding band structure. Specifi- 
cally, we calculate fundamental and collective modes and 
we compare the results with finite element analysis. In 
Section IV, we consider a specific example of nonlinear 
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FIG. 1: Antenna structure. 




dynamics of the system in which two resonant modes are 
coupled. This type of mode coupling, as we describe in 
detail, can allow detection of response in one mode by 
monitoring the coupled mode. 



II. DISCRETE MODEL 

A. Equations of motion 

The antenna-like structure we investigate consists of 
a central clamped-clamped diamond beam with dimen- 
sions L x W x h = 10.7/xm x 0.4^tm x 0.35/im and 40 
perpendicular cantilevers with dimensions I x w x h — 
0.5/im x 0.2/im x 0.35/im. The cantilevers are regularly 
spaced along the beam and symmetric with respect to the 
beam (see Fig. [TJ . There are N = 20 cantilevers on each 
side. In addition, the antenna is coated by a layer of gold 
whose length and width are those of the beam and whose 
thickness is t = 0.035/im. We denote by x the coordi- 
nate along the beam and by £ the transverse one. The 
deflection of the beam with respect to its clamps is called 
y(x, t) and the deflection of the jth lateral cantilever with 
respect to y(xj,t) is denoted by rjj(£,t). In the following 
we shall only consider vibrational modes symmetric with 
respect to the central beam. The deflection of the can- 
tilevers on either side of the beam are the same. They 
are described by the quantity r]j(£,t). Considering that 
all the elements of the structure are one-dimensional and 
using the Euler-Bernoulli beam theory (see for example 
Ref . [25 , 26 , 27 ) , the equations of motion of the structure 
have the form 



S t 



d 4 y 
'dx 4 ' 

d 4 i]j 
5^ 



m 6 S = -^E^ 



dt 2 
dt 2 



d£ 3 



S(x-Xj) (1) 



fh 



-fit 



d 2 Vj 
dt 2 ' 



j e {!,..., N} (2) 



where yj = y(xj, t) and together with the boundary con- 
ditions 

2,(0, t) = g(0, t) = y(L, t) = ^(L,t)=0 (3) 



»7i(0,t) 



(l,t)=0, (4) 



with j 6 {1, . . . , N}. The rigidities of the beam, £&, and 
the cantilevers, £t - that account for the presence of the 
gold layer - are respectively given by (see e.g. Ref. [26] ) 
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where Ed — 700 GPa and E g = 45.6 GPa are the mea- 
sured Young's moduli of diamond and gold, respectively. 
The masses per unit length of the beam and a cantilever 
are respectively given by 



= W{p d h + p g t) 



fit 



W 



fib 



(G) 



where p d = 4050 kg. m~ and p g — 19500 kg. m~ are the 
densities of diamond and gold, respectively. It is impor- 
tant to note that because of the smallness of E g with 
respect to E d , the rigidities of the beam and the can- 
tilever are barely affected by the presence of the gold layer 
(roughly 2.5 % higher). This is in contrast with their lin- 
ear mass that the gold layer increases by roughly 50 %. 
Thus, the presence of the gold layer merely increases the 
mass of the antenna without affecting its rigidity. 

The force density in the r.h.s of equation ([TJ represents 
the shear force densities exerted by the cantilevers on the 
beam. No momenta appear in this equation because, for 
modes that are symmetric with respect to the central 
beam, the momenta exerted by two opposite cantilevers 
cancel each other. The r.h.s. of equation ^ is due to 
the motion of the base of the cantilever that follows the 
motion of the central beam at £ = 0. 



B. Energy and Lagrangian 

The system of equations ([TJ and ^ (for j 
{!,..., N}) conserves the energy 
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(7) 



The first and second terms in Eq. ([T]) represent the elastic 
and kinetic energies of the central beam and of the 2N 
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cantilever, respectively. Similar to Hamiltonian ([7]), the 
Lagrangian of the antenna, Cd, is expressed as 



dx 



2 \dt 



dx 2 



N ' 
i=l n 



Ht ( drjj d yj 
2 I at dt 



£ t ( d 2 r) 3 



(8) 
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Eqs. Q-(|2j) can be established from the least action prin- 
ciple 5S d = 0, i.e. $g* = 
is given by 5 d = J £ d 



and 4^ = 0, where the action 



Vibrational modes 



The vibrational modes of Eqs. ([TJ-Q can be found as 
follows. Given the punctual nature of the force densities, 
Eq. reduces to 



dx 4 



d 2 y 
dt 2 



= for x e [0,L],x^ 
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(9) 



with j 6 {1, . . . , N}. At all points Xj, the function y(x,t) 
and its first and second derivatives are continuous. But 
its third derivative is discontinuous and its discontinuity 
is obtained by integrating Eq. (1) over an infinitesimally 
small interval centered around Xj. Thus, in addition to 
the boundary conditions given in Eq. ((3J, 4A other con- 
ditions apply: the continuity of y{x, t) and its first and 
second derivatives at x = Xj, j € {1, . . . , N}, and 
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= _2fi *^ (0 ** ) - (10) 

To simplify the system of Eqs. ([TJ-Q, as well as the 
boundary conditions, we shall work with the following 
non-dimensional quantities 
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Looking for mode solutions on the form y{x, t) — 
Y{u)cos{ujt) and r)j(£,t) = Hj(v) cos(wi), the equations 
to be solved become 



d A Y(u) 



du 4 



- ^ A Y{u) = for u € [0, l],u ^ Uj, (12) 



d*Hj(v) 



dv A 



1 4 (H j (v) + Y(u j ))=0, 



(13) 



and the boundary conditions Q and (10) read 

F(0) = F'(O) = F(l) = F'(l) = 0, 

Y(uj)=Y(uf) ; Y'(uj)=Y'(uf); 

Y"(uJ) = Y"(u+) ■ Y"'(uf) - Y"'(uj) = -flflf(0), 

H 3 (0) = i^.(0) = H'!{1) = H'!\l) = 0, (14) 



where Uj = j/(N+l), j € {!_ 



,N}. The solution to the 
system of equations ( 12 1-( |13[ ) is worked out in appendix 
A) From it, we eventually obtain the secular equation of 
the discrete model 



detM(w) = 0, 
where M(u) is a 2 x 2 matrix given by 



MM = T(/i) 



N 
3=1 



2aA 2 ( 7 )Kj( M )) 



(15) 
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In the expression above, 1 is the 4x4 unity matrix, 
a = w/W, A 2 (j) is defined in Eq. (30). The quantities 



[i and 7, defined in (111, are related to each other as 
fj, = "fL/l. They are the frequency dependent parameters 
of Eq. (16 1. Finally, the matrices Kj(^), TQu) and L are 
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; T( M ) = 


\ o 




V 



cos /i — sin fi\ 

sin fj, cos /i 

cosh fi sinh /i 

\sinh [i cosh /i / 



(17) 



where A T denotes the transpose of A and where Cj = 
cos(/iMj) ; Sj = sm(fiUj) ; chj = cosh(/ittj) ; shj = 
sinh(fmj) ; Uj — j/(N + 1). Using this, we obtain the 
exact solution of the discrete model. Nonetheless, the 
method is particularly tedious in practice and it does not 
allow for an easy analytical investigation of the problem 
unless the number of cantilevers to be treated is very 
small. The case TV = 1 is treated in appendix [5] and 
compared to the result given by the continuum approach 
developed hereafter. 



III. CONTINUUM APPROXIMATION 

A. Derivation of the model 

We now derive a continuum approximation for the 
system of Eqs. ([!])- Q. The idea is to "smooth out" 
the punctual force densities of Eq. ([!]) in such a 
way that the total shear force exerted by the jth can- 
tilever on the beam is the same but with a density 
that is continuous along the beam. The shear force 
exerted by the two cantilevers in the jth position is 
given by Fj — —2£ t (d 3 r]j/d£ l 3 )(0,t). Averaging its ef- 
fect over the interval [xj — S/2,Xj + 8/2] produces the 
force density fj — Fj/6 and filling the gap between 
the cantilevers requires 6 = L/N. Using this, we ob- 
tain a piecewise-constant force density given by fj = 
— (2N£ t / L)(d 3 j]j /d£ 3 )(0, t) whose continuous version is 
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straightforward if we now assume that the cantilevers 
form a continuum along the beam. Writing i]j(^,t) = 
rj(xj,£,t), the continuous version of the density simply 
reads f(x) = ~{2N£ t /L){d 3 i 1 /d£, 3 ){x, 0, t). Therefore, 
the equations of motion become 



d 4 y(x,t) d 2 y(x,t) 



dt 2 



c d 4 v (x,Z,t) d 2 V (x,Z,t) 

t-t fa 



2£tN d 3 V (x^,t) 
L d£ 3 

d 2 y(x,t) 



dt 2 



-IH- 



dt 2 



(18) 
(19) 



with clamped-clamped and cantilever-like boundary con- 
ditions for the central beam and the cantilevers, respec- 
tively 



y(0,t) = ^(0,t)=y(L,t) = ^ c (L,t)=0 



(20) 



dn d 2 v d 3 n 

V (x, 0, t) = q£& 0, t) = ^(x, I, t) =g^(x, I, t) =0. (21) 

This model is of course expected to provide better results 



as the density of cantilevers increases. Equations ( 18 I 
( 19 I conserve the following energy 



H r = dx 



2N 
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d 2 y 
dx 2 



fib 
2 



dx / d£ 



d 2 r\ 



2 



dr] 
di 



dy 
dt 



(22) 



Similar to the equations of the discrete model Q-([2]), 
Eqs. ( 18 1-( 19 1 can be deduced from the least action prin- 
ciple 5S C — 0, where iS c = JC C dt and where the La- 



grangian C c is obtained from Hamiltonian (22 1 by chang- 
ing the sign of the elastic energy terms. 

The continuum model presented above becomes more 
accurate as the density of cantilevers, N/L, increases. To 
get a sense of its limitations though, let us imagine that, 
for a given mode, the shape of the central beam has n — 1 
nodes, and consequently n "arches" . We then expect that 
the model holds if at least one cantilever per arch subsists, 
which basically yields the following condition of validity 
of the continuum model, 



n < N. 



(23) 



Indeed, should condition (|23j) fail to be satisfied, the "can- 
tilever continuum" would create some incrtial effects on 
an arch where no physical cantilever is to be found, which 
is not desirable. The limit of a single cantilever per arch, 
instead of two or more, can still appear rather arbitrary. 
But a comparison of the lower part of the spectrum of 
the continuum and discrete models for N = 1 done in 
appendix [B] reveals that their first frequencies are very 
similar. Therefore, thinking of the central beam as a 
collection of sub-beams (arches) pinned at the level of its 
nodes, it seems reasonable to take one cantilever per arch 
as the limit of validity of the continuum model. 



B. Vibrational modes 

Working again with the non-dimensional quantities de- 
fined in Eq. (11 1, we look for mode solutions of the form 



y(x, t) = Y(u) cos(wt) ; rj(x,£, t) = H(u,v) cos(wi). 

(24) 



The equations of motion (18) and ([19| become 
d 4 Y(u) 



du A 
d i H(u,v) 



d 3 H 
RN—„ - (u,Q) 



7 4 (H(u,v)+Y(u)) =0 



dv 4 

with boundary conditions 

no) = f(o) = y(i) = f(i)=o 



d 2 H 

dv 2 



(u,l) 



d 3 H 

dv 3 



(25) 
(26) 

(27) 

(u,l)=0. (28) 



Now, solving Eq. (26 1 yields 



H(u,v) = [^1(7) cos(7«) + A 2 (7) sin(7w)+ 

^3(7) cosh(7i;) + ^4(7) sinh(7w) — 1] Y(u) 
= H(v)Y(u), (29) 

where the coefficients ^(7), determined by the boundary 
conditions ( 28 1 , are given by 



^i(7) 
^2(7) 
^3(7) 

Mi) 



1 + cos 7 cosh 7 — sin 7 sinh 7 

2(1 + cos 7 cosh 7) 
cos 7 sinh 7 + sin 7 cosh 7 

2(1 + cos 7 cosh 7) 
1 + cos 7 cosh 7 + sin 7 sinh 7 
2(1 + cos 7 cosh 7) 
cos 7 sinh 7 + sin 7 cosh 7 
2(1 + cos 7 cosh 7) 



(30) 



From these results we can calculate the r.h.s. of Eq. ( 25 ) , 



-RNd 3 H/dv 3 (u,0) = 2RN"/ 3 A 2 (j)Y(u), and we finally 
obtain 



d 4 Y{u) 
du A 



P a Y(u) = 0, 



(31) 



with 



/ 3 4 = / + 7 3 ^ COS7Sinh7 + Sin 7 C ° Sh7 . (32) 
1 + cos 7 cosh 7 



Eqs. (31 1 and A32h show that our continuum model re- 



duces to a simple clamped-clamped beam equation with 
a parameter /? that depends in a non trivial way upon 
the mode frequency lu. Indeed, according to Eq. (11), 



the parameters /i 4 and 7 4 are quadratic functions of the 
frequency and the parameters R and N are simply con- 
stants. Now, because Eq. (31 1 along with boundary con- 



ditions (27 1 is merely a clamped-clamped beam equation, 
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it can be solved to obtain 

Y(u) = A{ cos(fti) — cosh(/3u) 

cos /3- cosh P\ . \ 

(sm(pu) — smh(pu)) >, (33) 



sin /3 — sinh (3 
with the secular equation 

cos(ft cosh(ft = 1. 



(34) 



Frequency spectrum of the continuum model: 
Band structure 



We see that the shape of the vibrational modes of the 
central beam is uniquely determined by the parameter (3 
solution to Eq. (34 1. The solutions ft can be evaluated 
numerically: ft = 4.730040745, ft = 7.853204624, ft = 
10.99560784, and ft ~ (n+l/2)7r as n is large. To a given 



ft corresponds a single modal shape for the central beam, 
Y n (u), but an infinite number of modal shapes i/„ fc(u, v) 
and frequencies w n ,fe, k > 1, obtained by solving Eq. (32 1 
for w . To see that the number of solutions to Eq. ( 32 1 
for P = ft is infinite, we first rewrite it in terms of the 
variable 7 by using fi = jL/l and denoting by = 
L/ib and = //it the masses of the beam and of one 
cantilever, respectively, we finally find 



Q(7.A0 



1 



cos 7 



1 



^ft 
L 7 



valid for cos 7 cosh 7 
when 7 is large, Eq. (35 1 reduces to cos 7 



cosh 7 

2Nmt cos 7 tanh 7 + sin 7 

m b 7 

1^0. 



(35) 



It becomes now clear that 
0. This equa- 
tion has an infinite number of solutions with asymptotic 
behavior j n> k = (k — l/2)w as k — > 00. Interestingly, these 
asymptotic solutions are independent of all the parame- 
ters entering Eq. ( [35] ) , and in particular, of ft provided 
the latter is finite, i.e. n 00. Finally, the frequency 
corresponding to a given value j„ t k is given by Eq. (Ill 
as 



(36) 



To understand the structure of the spectrum, we have 
plotted in Fig. [2] the frequencies uin.k normalized to the 
fundamental, u>n, for N = 20 and the antenna dimen- 
sions indicated in |II A| As we can see, the spectrum con- 
sists of "bands" separated by gaps. The band number 
is the label k that we have attributed to the solution of 
Eq. (35 1, 7 rii fc. To elucidate the appearance of bands, we 



note that as n increases, the solution to Eq. (34 1 becomes 
large and is asymptotically given by ft ~ (n+l/2)n. For 
a finite solution, 700. fe, of Eq. (32 1 to exist as n — > 00, 



the denominator of the second term of the r.h.s of (32 1 
needs to vanish. Hence, 
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FIG. 2: Frequency spectrum obtained from Eq. (35 1 for 
iV = 20 and the antenna parameters given in |II A| Solid lines 
correspond to n < N (physically relevant frequencies) and 
dashed ones to n > N. 



This equation is the well-known secular equation for a 
simple cantilever. It provides here the upper band edge 
of the kth band whose frequency, determined by Eq. ( 36 ) , 
reads 



^00, k 




(38) 



~ 100.103, 

l/2W 7 i,i) 2 



1 + cos(7oo,fe) cosh(7oo,fc) = 0. 



(37) 



For the case at hand, we have o; 00i 
^00,2/^1,1 — 627.336 and LJoo : k/^i,i — ((k - 
with 71.1 ~ 0.187413 for k > 3. It is easy to see that the 
accuracy of the latter formula improves exponentially as 
k increases. For k — 3 the relative error to the exact 
result is already as small as 0.02%. An approximate ana- 
lytical expression for 7^1 is given in the next section while 
results for the lower edge of the fcth band, 71 k, k > 2, 
will be given in section |III E| 

At the edges of each band but the lower edge of the first 
one that we shall refer to as the "fundamental band" , fre- 
quencies clearly accumulate. However, close to the lower 
band edge, a finite number of frequencies cluster while 
an infinite number accumulate at the upper band edge. 
Each band contains all possible values of ft, 1 < n < 00, 
i.e. all possible modal shapes for the central beam, Y n (u). 
Even though modes with the same shape Y n (u) repeat- 
edly appear within each band, they differ from band to 
band because their cantilever continuum, H„ t k(u,v), de- 
pends on 7„ fc that is, on both n and k (see Eq. (29)). 
Note also that, within a given band, the frequency in- 
creases with the excitation level of the central beam. 

At mid-band, modal frequencies are somewhat sparse 
and more regularly spaced than at the boundaries. This 
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FIG. 3: Frequency spacing, (w n +i,fc — uo n ,k) /ui\.i, versus fre- 
quency (in red) for the same parameters as in Fig. [5] In blue, 
analytical results obtained from Eq. (39 1. Curves are solid 
for n < TV and dotted for n > TV. 



is confirmed in Fig. [3] that displays the normalized spac- 
ing between two consecutive frequencies versus the nor- 
malized frequency itself. Notice that, the inverse of this 
function is nothing but the normalized "density of states" 
of the spectrum. Mid-band frequencies can be approxi- 
mately evaluated once we note that the function ^2(7) 
appearing in secular equation ( 32 1 is small away from 
the boundaries 700, & given by Eq. (37) and increases 
slowly and regularly in this region. Rewriting Eq. (32) as 
fi 4 — f]^ + 4aNA 2 ((j,l/L)fj? = 0, we seek a solution close 
to ji *~ (3 n . With a first Newton iteration we obtain 



In.k ^ 7. 



(0) 



aNA 2 { 1 ( ° ) ) 



(3 n + 3aNA 2 ( 1 i;") + aN(3 n j;A' 2 ( 7 k 0) ) 



/4I « 



LA k 
21 



(39) 



where 7^ = a = w/W and A 2 (~f) is given in Eq. 
pfjj ). In the last equation, A fc = 700,^ — 71^, is the 
width of band k and /it is determined as the solution to 
A 2 (n* k l/L) = 0. For k large enough, /xjj; ~ (k - l/4)nL/l, 
k > 1. In this case, we have also (3 n ~ (n + l/2)n, n > 1. 
Approximating the bandwidth by ~ 700, & — 7oo,fc-i — 
7r as k is large enough, the condition of Eq. ( p9| ) becomes 
\{n+ 1/2)1/ L- (k- 1/4)| < 1/2, which provides bounds 
on the beam excitation (n) of band k for which (39 1 is 
valid. In this case Eq. (39 1 simplifies to 



p n - aNA 2 



\Pn-(4\ < 



21 



(40) 



Although simple, this last expression is not very accurate 
for the first bands. 

Very few of the many frequencies occurring in the spec- 
trum of the continuum model can be detected experi- 
mentally. This is mainly due to the modal response of 
the structure which, for a given driving power, is drasti- 
cally suppressed as its excitation level increases as we 
shall prove in the study of a damped driven antenna 



done in section IIIH Moreover, the antenna is generally 
driven by a harmonic force that is constant over the beam 
which does not allow for the observation of antisymmet- 
ric modes, Y n (u) with n even. In what follows, we show 
how to evaluate the frequency of the fundamental mode 
of the structure and explain the clustering phenomenon 
observed at both edges of each band. 



D. 



Small 7 solution, fundamental mode and mass 
loading 



If a small solution to Eq. (35) exists, 7 < 1, we can 



obtain its approximate expression by expanding Eq. ( 35 1 
around 7 = 0. In this case we find that 



7 



1/4 



IJ3 

T' 



(41) 



where fi a — (m^ + 2Nm t )/L is the mass per unit length 
of the entire antenna (beam plus cantilevers). Replac- 
ing this result in Eq. ( [36]) and using the fact that 
(Ht£b) I \Hb£t) — 1 ( see Eqs. (5 1 and ((6|)) yields 



(42) 



In this limit, the inertia of the cantilevers is negligible and 
the frequency of the vibrational mode is the frequency 
of a clamped-clamped beam whose mass per unit length 
includes the mass of the central beam plus the mass of the 
cantilevers. According to Eqs. (29 1 and (30 1, the mode 




shape H (u, v) corresponding to this small 7 solution is 
found to be 



H(u,v) 



7 

— 1 
24 



l (v 2 -4v + 6)Y(u). 



(43) 



Notice that this function is proportional to 7 4 and is very 
small. This confirms that the cantilevers hardly move at 
all and that their mass only loads the central beam. Fi- 
nally, given our parameters, the condition 7 <C 1 is typi- 
cally satisfied for the first positive root of Eq. ( 35 1 when 
j9 = P\. It is then valid for 71, 1, i.e. for the fundamental 
mode. Indeed, using Eq. (41 1 and the antenna param- 

- 0.1874150... while 
1 = 0.1874136..., i.e. 



we find 71,1 



eters of section II A 
the exact numerical value reads 71 
a relative error of 7 x 10 -4 %. 



E. Frequency clustering 



The frequency clustering evoked in |III C[ which has also 
been observed numerically in 3D finite- clement Simula- 
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tions, can be explained with the help of secular equa- 
tion ( 32 1 . To evaluate the accumulation of frequencies 
occurring at the upper edge of the kth band (c^oo.fc) 
and at the lower edge of the (k + l)th band, we seek 
a perturbative solution of Eq. J32j) around 7 ~ 700, fc 
given in Eq. (37 1. Using Eq 

<2(7oo,fc, Pn) / Q' (loo, k, Pn) 

k+ 1, and Q' is the derivative of Q with respect to 7, we 
obtain 



35]) and i nk ~ 700, k 
where k stands cither for k or 



7„E - 7oo,fc + ' 



OQ.k 



2aNl \t 



7oo,fc 



-th^ 
-th~, 



-th~, 



th„ 



1- 
1 

7oo,fc 



\^7oo,fc 

1 



(44) 



where i~ = tan 7 and th-y = tanh7. This solution is valid 
as far as the corrective term is very small compared to 
7oo,fc- This leads to the two conditions below. 



1. Upper band edge 

With respect to the frequency clustering at the upper 
band edge of the kth. band, Eq. (44 1 starts to be valid as 



soon as p n is large enough. In this case, it is easy to see 
that this approximately leads to the condition 



Pn > 



£70 



1 



2aNl\ 



1/4 



L 



(45) 



This condition can always be satisfied for large enough n, 
given that in this case, p n ~ (n+l/2)ir (for simplicity, we 



have discarded the term (t^ 



■*' l 7 c*)/( fc roo,*+*' l 7oo I *) 



that is close to unity). The longer the cantilever, the 
easier the condition (451 to be satisfied. For an antenna 



structure with high cantilever to beam length ratio, we 
may expect to observe a clustering of modes around the 

Experimentally, however, 



frequency determined in (38 



this very much depends on t 



le sensitivity of the measur- 



ing device given that highly excited modes for the beam 
are difficult to detect. Moreover, the range of validity of 
our model, determined by the condition n < N, restricts 
the possibility of observing the beginning of an upper 
edge clustering to the fundamental band. 

Provided (45 1 is satisfied, Eq. (44 1 simplifies to 



ln,k ~ 7o 



Pi 



th~. 



Pi 



.k + ^Too,* 

. + *T 

00 , k Joe 



(46) 



The second equality has been obtained from ( 37 1 . As we 
can see, the corrective term to 700, & is always negative 
which confirms that the frequencies are indeed accumu- 
lating at the upper edge of the band as n increases. It 
is also possible to show that the (adequately normalized) 
modal shape of the cantilever corresponding to Lo n ^ fac- 
torizes asymptotically as the product of the feth mode of a 



simple cantilever by the nth mode of a clamped-clamped 
beam, H n ^ k (u,v) ~ 4>k{v)Y n (u) as n — > 00. Thus, as the 
excitation of the beam becomes high enough, the modal 
shape of the antenna is given by a simple cantilever mode 
modulated by the profile of the central beam. The same 
result is obtained for a weakly dissipative driven antenna 



2. Lower band edge 

A frequency clustering at the lower edge of each band 
does not always occur. This essentially depends on the 
geometry of the antenna. For the correction given in Eq. 
di4l to be valid, 



Pn < 



I 



2aNl ( th laa h + t la 



L 



th„ 



-1 1/4 



(47) 



has to be satisfied. For large cantilever to beam length ra- 
tio, however, inequality (47 1 cannot be satisfied because 



the square-bracketed term becomes negative. This means 
either that expanding the secular equation around the up- 
per edge of the A:th band does not provide any reliable 
information on the lower edge of the (k + l)th band or 
simply that no clustering occurs at the lower band edge. 
For the geometry of the antenna described in the intro- 
duction, frequency clustering occurs and the band gap 
^■k,k+i between the kth and (k + l)th bands can be ac- 
curately evaluated with the help of Eq. (44 1 as 



A fe . 



fc+i 



7i,fc+i — 7oo,fe- 



(48) 



This expression is accurate to less than 5% for the first 
gap and to less than 0.5% for the second one and its 
accuracy improves drastically as k increases. 



F. Energy and Lagrangian 

1. Modal energy 



Using expressions (p} and ( |24p4| ), we can evaluate the 
energy E n<k of the mode ip njk = {y n (x 1 t),r] n ^(x,^t)) T . 
After some algebra, we find 



E n ,k — 



Mn,kA'<k 



where 



M n , k = m b + 2Nm t L n u. 



(49) 



(50) 



The quantity L n ^ k , that takes into account the dynamics 
of the cantilevers at the level of the effective mass M nyk , is 
given by L n<k = L(j n>k ) where, defining h(v) = H(v) + 1, 



1 ( cos 7 + cosh 7 \ 2 3^2(7) 



LOy) = / h 2 (v)dv 
KU J Q w 4 VI +COS7COSI17/ 2 7 



(51) 
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Incidentally, we also note a relation that will prove useful 
when dealing with the motion of a driven weakly damped 
structure in section ITlI H 51 



i d 

2^3 



(52) 



Eq. (49 1 is exactly the energy of an effective harmonic 
oscillator with amplitude A, frequency u>n : k and mass 
M„.fc. Notice that, according to our definition of Y(u) 

in (33 1, the amplitude is defined as A 2 = J^Y 2 (u) du 
which is independent of n and k. Other definitions of 
the amplitude (like the maximal deflection of the cen- 
tral beam at the mid point, for example) are typically 
mode-dependent and redefine the effective mass of the 
structure. As we see, the effective mass of the antenna 
is renormalized by the factor L n j. that affects the total 
mass of the cantilevers, 2Nm t . In the specific case of the 
fundamental mode, we have shown in the previous sec- 
tion that 71 1 <C 1. Therefore, L\ \ ~ 1 and the effective 
mass of the antenna, M\ t \ ~ m b + 2Nmt, is its actual 
mass. This confirms the finding of section III D that in 



the fundamental mode the cantilevers have a pure mass 
loading effect. We can also use Eq. (49) to define an 
effective mode-dependent spring constant K n ^ as 



(53) 



Once again, this result depends on the definition chosen 
for the amplitude and is valid here for A 2 = f^Y 2 (u) du. 



'nk 
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FIG. 4: Ratio, r„ t k (see Eq. ( 56 1 ) , of the maximal elastic 
energy of the beam to the total energy of the antenna versus 
the normalized frequency, uj n ,k/ui,i- Antenna parameters are 

20.' 



given in section [II A| and TV = 20. Solid lines indicate physi- 
cally relevant results for the continuum model (n < TV) and 
dotted lines the results for (n > TV). 



2. Energy partition 

To complete the modal analysis of the antenna, we 
evaluate the energy partition of each mode that is, we 
compare the energies of the cantilever continuum and of 
the central beam to understand which part of the struc- 
ture is the most active. For a mode, the elastic energies 
of the central beam and the cantilever continuum reach 
their maximum at the same time and at this point their 
kinetic energy is zero. The sum of the maximal elastic 
energies of the beam and the continuum is then equal to 
the total energy of the system. We thus choose to calcu- 
late the ratio, r n ^, of the maximal elastic energy of the 
beam, U b] n,ki to the total energy of the structure, E n k , 
as a function of the frequency of the mode, Lo llt k- The 
function r n ^ ranges from to 1 and is expected to be 
useful in the design of antenna structures with a specific 
task. The central beam elastic energy is given by 



St fd^y 
2 I dx 2 



U h [y{x,t)] = J dx 
o 

Using expressions ( 24p4 ), we obtain 



(54) 



Ut-n^k = ma,xU b [Y n (x/L)cos(u; n ^t)} = —. ^A 2 (3^ 



Nm t A 2 tu 2 



"* \2Nm t 



2U> 

, A 2 (jn,k 
In.k 



(55) 



Using expression (49 1 for the total energy E n ^, we finally 
have 



Ub-n,k 
E n k 



2Nm t 



2Nm t 



f L n .k 



(56) 



where L n ^ is given in Eq. (51 1. In Fig. [4] we display r n ^ 
versus the normalized frequency oj n ^/^i,i- We observe 
that at the edges of all the bands but the lower edge of 
the first one, r n ^ almost vanishes. This means that most 
of the energy of the structure is located in the cantilever 
continuum. By contrast, in the middle of the bands, most 
of the energy is stored in the central beam. This confirms 
the results of the previous sections showing that close 
to the band edges, the parameter 7 almost satisfies the 
"cantilever" secular equation 1 + cos 7 cosh 7 = 0, which 
clearly indicates that the spectrum is mainly governed by 
the dynamics of the cantilever continuum in this region 
while in the "middle" of the bands, j n .k ~ (l/L)(3 n (see 
Eq. (|40|), which makes the spectrum close to that of a 
clamped-clamped beam. The peculiarity of the first band 
lies in the fact that for j n> i close enough to zero, r n> \ ~ 
1-7^ 1 /[20(a-+l)] where a = m b /2Nm t . The energy of 
the very first modes and, in particular, the fundamental, 
is thus essentially stored in the beam. This confirms the 
results obtained in section IIII Dl 



10 



3. Lagrangian 



/i t 2 — 1-55 GHz, within 3 percent of the simulated value. 



It is also interesting to evaluate the Lagrangian C n ,k{t) 
of the system and in particular its time average that we 
will use in the next section to evaluate the effects of small 
nonlinearitics and damping on the structure. We find 



M n k A 2 ujl k 

— cos(2w„ !fc t) =*> (£ n<k (t))=0, 

(57) 

where (...) denotes the time average. The modal La- 
grangian of the linear continuum model is then similar 
to a harmonic oscillator's and its time average is zero. 
This result will be used in the perturbative treatment of 
a weakly nonlinear and dissipative antenna in section |IV| 



G. Fundamental and first collective modes: 
comparison with finite-element results 

In this section, we compare the results of our contin- 
uum model (CM) to those obtained with a finite element 
method (FE) that treats the vibrations of the antenna 
within the frame of three-dimensional elasticity theory. 
We are primarily interested in the fundamental and first 
collective modes, that is in the first modes of the first 
and second band, respectively. These modes are easy 
to observe experimentally. They are related to the pa- 



rameters 7i.i and 71.2, respectively. Using Eq. (35 1 and 
the antenna parameters given in |II A we find that the 
two first roots of Q("f, with /?i = 4.73004... are 
71,1 = 0.187413 ... and 71,2 = 2.046440 .... If we use 
Eq. pll, one finds 71,1 ~ 0.187415. It shows that this 



approximate solution is very reliable for the fundamental 
frequency. The corresponding frequency, obtained from 
Eq. ([36f, is / M = w m /2tt ~ 24.7 MHz, slightly higher 



than the frequency observed in simulations (23.6 MHz). 
This might be explained by the effective cross-sectional 
stiffness that we approximate for the two material lay- 
ers. The frequency of the first collective mode calculated 
from 71,2 = 2.046440 is / 1)2 ~ 2.94 GHz which is much 
higher than fff — 1.51 GHz, the frequency from finite el- 
ement simulations. A reason for this discrepancy is that 
all elements of the antenna structure are considered as 
one-dimensional in our model. Consequently, their length 
is the only dimension taken into account in the dynam- 
ics of the system. However, in the real structure, the 
central beam has a nonzero width, W, on the order of 
the cantilever length, I. For motions of the cantilevers 
comparable to the beam's, the shear momentum they ex- 
ert on both of its sides becomes large enough to bend it 
laterally with respect to its mid-line, an occurrence in- 
deed confirmed by our three-dimensional finite-clement 
simulations. In first approximation, we can account for 
this effect by assigning to the cantilevers an effective dy- 
namical length ranging from their actual length, I, for 
small amplitude motions to I + W/2, for large ones. If we 
carry out the substitution / — > / + W/2 in Q(j,(3i) and 
reevaluate 71,2, we find 71,2 = 2.085115 and a frequency 



The shapes of the central beam predicted by the con- 
tinuum model are the same for the fundamental and first 
collective modes. They are given by Y\(u), obtained from 
Eq. ( 33 1 for (3 = j3\. The shapes of the cantilever con- 



tinua, however, are different and given by equation (29 1 
as Hi t i(u,v) and Hi^{u,v), respectively. A comparison 
of these results to those of the finite element simulation 
is done in Fig. [5] The FE and CM results for the funda- 
mental and first collective modes are given in the upper 
left and upper right panel, respectively. As indicated, 
the FE result is to the right of the CM result. In both 
cases, the color code (online) indicates the displacement 
(amplitude) of the elements with respect to the clamps 
of the central beam (i.e. its extremities). For the CM 
results, only half of the absolute displacement of the can- 
tilever continuum, Yx(u) + Hi t k(u,v), k = 1,2, v > 0, is 
displayed. For v — 0, this displacement is precisely the 
beam's, Y\(u). As we see, the continuum model predic- 
tions are in excellent agreement with the finite element 
results. For both the fundamental and the first collective 
modes, the shape of the central beam is in the funda- 
mental mode of a bare beam, Yi(it). In the fundamental 
mode, the deflection of the cantilever continuum with re- 
spect to the beam is imperceptible as predicted by the 
continuum model in the small 7 limit (see HID ). Indeed, 
from Eq. |43|, we have Hi t \{u, v) cx 7i 1^1 (u), which 
means that the continuum deflection is roughly 5000 
times smaller than that of the beam. This is markedly 
different for the first collective mode where the motion of 
the continuum is on the order of the beam amplitude as 
seen on the upper right panel. Note that, if their ampli- 
tudes are very different, the actual shapes of the continua 
for the fundamental and first collective modes are similar 
and are given to a very good approximation by the fun- 
damental mode of a bare cantilever. This somewhat sur- 
prising fact becomes clear on the modal expansion of the 
cantilever continuum, H nt k(u,v) = H n _k(v)Y n (u). Intro- 
ducing the normalized cantilever modes, ipi(v), satisfying 



Too i^ii v )^ we obtain from Eq. (26 1 



H„. k (v) 



H 7 4 

1=1 ,oc M 



(58) 



where ^1 = f Q ipi(v)dv = — (2/7) (cos 7 + cosh 7) /(sin 7 + 
sinh7), where 7 = 7oo,;- Now, we have seen in section 
HIE that, for k > 2, is very close to 7oo,fc-i- Accord- 
ing to Eq. ( 58 1 , this means that the mode selected for 
the cantilever continuum is ipk-i(v). In particular, for the 
first collective mode, 71,25 the continuum adopts basically 
the shape of the fundamental cantilever mode, ip\{v). For 
the modes of the fundamental band (k = 1), this is also 
the fundamental cantilever mode that is predominantly 
excited, even if its amplitude is so small compared to the 
beam's that the cantilevers seem flat at the scale of the 
figure. 

To analyze further the similarities between the FE and 
CM results, the first excited modes of the fundamen- 
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FIG. 5: Comparison between finite element (FE) and continuum model (CM) results for the fundamental mode (upper left 
panel), the first excited mode of the fundamental group (lower left panel), the first collective mode (upper right panel), and the 
first excited collective mode (lower right panel) . The parameters for the antenna structure are those quoted in II A 



tal and collective bands are displayed in the lower left 
and right panels of Fig. [5] respectively. We clearly ob- 
serve in both cases that the beam deflection is in the 
first excited mode of a bare beam, Y 2 (u), as predicted by 
the continuum model. Moreover, as for the fundamen- 
tal mode, the cantilevers of the first excited mode of the 
fundamental group are hardly moving. Evaluating the 
first root of Q(>y, (3 2 ) = with /3 2 = 7.85320 . . . yields 
72.1 = 0.31114.... This is small enough for the results 
of section HID to hold. Indeed, expression (41 1 gives 
72,1 ~ 0.31116.... Therefore, the continuum deflection 



can safely be evaluated from Eq. (43) and we find that 



it is roughly 100 times smaller than the beam deflection, 
which explains that no cantilever motion can be detected 
in the FE results. Once again, this is markedly different 
for the first excited mode of the collective group whose 
cantilevers, according to the finite element results, are 
experiencing a deflection comparable to the beam ampli- 
tude in full agreement with the predictions of the contin- 
uum model. Note that for this two first excited modes 
of the fundamental and collective bands, the cantilever 
deflections adopt the shape of the fundamental mode of 
a bare cantilever, ipi(v), as indicated by the modal ex- 
pansion (58 1. 



Even though one-dimensional in essence, the contin- 
uum model gives a good qualitative understanding of the 
modes of the 3D antenna. It correctly reproduces the 
modal shapes of the central beam and the cantilevers ob- 
served in the finite clement simulation and is able to ex- 
plain the frequency clustering occurring in the spectrum 
of the structure. We use it in the next section as the 
basic model in the investigation of the effect of a two- 



frequency driving on the response of a weakly nonlinear 
and dissipative antenna. 



H. Driven damped system: Exact solution 

We conclude this section devoted to the linear system 
by calculating the exact solution of a damped antenna 
structure driven by a spatially uniform harmonic force 
density, f(t), in the continuum approximation. This force 
density can, for instance, be exerted by an external mag- 
netic field, B, orthogonal to the flexural vibrations of the 
central beam, acting on an ac-current with frequency ujd 
passed through the thin layer of gold coating the struc- 
ture. The beam vibrations then generate in turn an elec- 
tromotive voltage V em f(t) at the clamped ends of the 
antenna that is proportional to the rate of change of the 
magnetic flux, 



V emf (t) = B J 



dy(x,t) 

at 



dx. 



(59) 



Ultimately, this voltage can be monitored to determine 
the motion of the central beam of the structure and 
to derive its "spectrum" i.e., the time root-mean-square 

of the induced voltage, (V^ rn j.-(t)) , versus the driving 

frequency, u>&. This detection scheme has been used 
in the previously reported experimental measurements 
of the antenna resonators in Ref. P2J Elj- Other ac- 
tivation/detection schemes such as electrostatic, piezo- 
electric, and optical are typically also sensitive to the av- 
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erage transverse displacement of resonating elements, and 
the following analysis applies quite generally. 



Model 



Supplementing the equations of motion (18)-(19l with 
damping and driving terms we obtain 



£, 



£, 



d A y d 2 y 



dx 4 
d 4 i] 



Mb 



dt 2 
d 2 r] 



dy 
dt 

drj 
di ' 



2£ t N d 3 r) 
L W 
d 2 y 



dt 2 



dy 
dt 



/(*) (60) 



(61) 



Boundary conditions for the beam and cantilever deflec- 
tions are the same as in (20 1 and (21 1. The harmonic 
force density is given by fjt) = f cos(w<ji)- Note that 
in Eq. (611, the damping term affecting the cantilevers 
involves their absolute displacement, y{x,t) + rj(x,^,t), 
rather than their relative displacement, r](x,(;,t). The 
choice of an appropriate damping term depends of course 
on the type of damping experienced by the structure. To 
simplify, we consider here that damping occurs through 
air friction and is then proportional to the absolute ve- 
locity of the cantilevers. Moreover, it is proportional to 
the surface in contact with the ambient air. For that rea- 
son, the damping per unit length is proportional to the 
width of the element involved and then vt/^b — w/W = 
Ihl^b — £t/£b- Material damping, whether of the vis- 
coelastic or hysteretic type (see for instance Ref. [2"4"]). 
would essentially affect the rigidities of Eqs. (60 1 and 
(61 1 in such a way that £ is replaced by £ + £*dfdt. The 



exact solution of system ( 60 )-( 61 1 can still be obtained in 
this case. 



2. Exact solution 

Introducing the Fourier transform g(u>) = J dte luJt g(t) 
and its inverse g(t) = J dui e~ luJt g(ui) /(2ir) and using the 
dimensionless quantities defined in (fTTI) , we can cast Eqs. 



(|60J-(|61| into 

d*y 
du 4 ' 
d 4 fj 
dv 4 ' 



4 ~ 



-NR 



d 3 rj 
dv 3 



v=0 



4 ~ 4 _ 



(62) 
(63) 



where F(t) = Fcos(iVdt) with F = L 4 f/£b and 
/ 4 



7 C 



Si 
S b 



((i t Lo 2 + iv t u), 
(fj,bW 2 + ivyj) 



Lie 
I 



(64) 



These last parameters are the complex version (because 
of the presence of damping) of the parameters \i and 7 
defined in (11). Note also that F has the dimension of a 



length. As in the earlier case, we can solve Eq. (63 1 with 
its boundary conditions: 



fj(u,v,u) = H(v,uj)y{u,uj), 



(65) 



where 



(66) 



H(v,u) = [yli(7 c )cos(7 c u) + A 2 (7 c )sin(7 c i;) 
+A 3 (7c) cosh(7 c w) + A 4 (-f c ) sinh(7 c u) - 1] 

The coefficients Ai are the same as those given in pO] ). 
Reinstating this expression in Eq. (|62[), we finally obtain 



du 4 



p 4 c y(u,u) = F(u;), 



(67) 



-)' 3 RN C0S lc 8inh 7c + Sin 7c C ° sh 7c (68) 
1 + cos 7 C cosh 7 C 



where 



which is the complex analog of the secular equation (|32j) . 
Finally, applying the appropriate boundary conditions to 
y(u, ui) we can, after some algebra, cast the solution to 
Eq. m\ into 



y(u,uj)= T(/3 C ) 



cos(/3 c (n - 1/2)) 



- 1 



sin(/3 c /2) 

cosh(/3 c (u- 1/2)) 
+ sinh(/3 c /2) 

where we have defined 

tan(/3 c /2) tanh(/3 c /2) 
tan(/3 c /2) + tanh(/3 c /2) 



/3 C 4 ^ 



(69) 



(70) 



Now, the driving f(t) being harmonic, we have F(lo) — 
Ftt[5(lj— uj^+diuj+uj^)]. Noticing that the sign inversion 
u> — > —lu amounts to taking the complex conjugate, we 
can finally show: 



y(x,t) = FRe{exp(-iud,t)y h (u,Ud)}, 
where yh{u,uj) — y{u,uj) / F{lo). Then 



(71) 



y(x,t) = FRc 



-iudt 



in 



T{(3 C 



( cos(/?c(m 



1/2)) 



cosh(/3 c (w - 



sin(/3 c /2) 
1/2)) \ 



sinh(/3 c /2) ) 



(72) 



Expression ( 69 1 makes it clear that the shape of the cen- 



tral beam induced by the force density f(t) is symmetric 
with respect to its midpoint, u = 1/2. This is expected as 
the force density itself possesses this symmetry. Conse- 
quently, none of the antisymmetric modes of the central 
beam are excited by this method. Moreover, in pres- 
ence of dissipation, the beam shape never corresponds 
to an exact symmetric modal shape even when the driv- 
ing frequency is one of the structure modal frequencies. 
We can see that for weak dissipation, however, the de- 
nominator of T(/3 C ) becomes small (on the order of v t ) 
when LUd ~ ^>in+\,k- This is because the solutions to 
tan(x/2) + tanh(x/2) = are precisely the /?2n+i's. This 
eventually leads y(x,t) to assume a shape close to the 
mode Y 2n +i(u). 
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3. Modal expansion 

As we have seen earlier, the modes of the antenna 
structure are such that the central beam possesses the 
exact shape of a clamped-clamped mode, Y n (u) (see Eq. 
p3|)). However, the force f(t) applied to the beam ex- 
cites now all symmetric clamped-clamped modes. To get 
a sense of which modes are predominantly excited, we 
first expand the Fourier transform of the deflection as 
y(u,uj) = *5Z n y n {w)Y n (u), insert it in Eq. (67! and solve 
for y n (tjj). The beam deflection given by Eq. (71 1 then 
reads 



y(x,t) = F 



OO 

E 

n=0 



r 2n+ i COS (u> d t + 02n+l(<4i)) 



£ 4 



Y 



2n+l 



2n+l | 



(73) 

where T n = J Q Y n {u)du and where $ n (w) = arg(/3^(w) — 
/3*). The reason why the sum runs exclusively over odd 
numbers in Eq. ( [73] ) is because integrals of antisymmetric 
clamped-clamped modes, I^n, vanish which confirms that 
none of them is excited by f{f). For symmetric modes, a 
simple calculation yields r2„+i = 4tan(/?2n+i/2)//32n+i- 
Clearly, expression ( p73~| ) shows that the mode Y2n+l( u ) is 
singled out when the driving frequency is close to one of 
the modal frequencies (J2n+t,k an d the dissipati on is we ak 
enough (\(3^(ujd) — Pzn+il K v t m this case, sec III H 5 1. 

Also of interest is the modal expansion of the cantilever 
continuum. Introducing the normalized cantilever modes, 

a tk ipk(v), we can 

and finally solve 



ipk(v), k > 1, satisfying i>% v \v) = 7; 



expand H(v,oj) in (66) as H(v,oj) — 
uses Eq. ( 65 1 and insert it in Eq. ( 63 1 



for H k (uj). This yields 



OO 

E 

k=l 



■7 C 4 H 



7 4 t 

' OO, ft! 



V , fe(^) I 



(74) 



where = L ipk(v)dv = — (2/7)(cos7+cosh7)/(sin7+ 
sinh7), where 7 = 7^^. Hence, ^ fc ~ -2/7^, as fc be- 
comes large. The interesting point about this calculation 
is to show what cantilever mode is selected according to 
the driving frequency. When the latter is close to the 
modal frequency w n k with n small, and provided the dis- 
sipation is weak enough, 7 c (w„j c ) ~ 7n,fe- Now, we have 
IIIEI that, for k > 2, 7^ is very close 
this means that for 



(741 



seen m section 
to 7oo,fc-i- According to Eq. 
tUd ~ the mode selected for the cantilever contin- 

uum is ipk-i(v). In particular, if the system is driven 
near the first collective mode frequency, UJ1.2, the contin- 
uum adopts basically the shape of the fundamental can- 
tilever mode. Note that, for driving frequencies in the 
fundamental band (fc = 1), this is always the fundamen- 
tal cantilever mode that is predominantly excited. For 
other bands (fc > 2), the continuum interpolates between 
the shape of the mode ipk-i( v ) close to the lower band 
edge and ipk{v) close to the upper band edge. In prac- 
tice, however, as our model provides consistent results for 
n < N only, the range of physically relevant frequencies is 
restricted to the lower band edge. For all peaks but those 



of the fundamental band, therefore, the relevant modes 
of the structure at frequency u> n ,k are Y n (u) and ipk-i(v)- 

4- Amplitude- frequency spectrum 

Most of the time, this is not the beam deflection y(x, t) 
but rather its average over the beam length (or its time- 
derivative, see below) that is detected. For that reason, 
we define 



y{x, f)dx. 



Using expression ( |72| , we immediately obtain 
y(t) = A(uj d ) cos(uj d t - 0{uj d )), 



(75) 



(76) 



where the amplitude A(u>) and the phase 9(u>) are respec- 
tively given by 



9(u) 



F 

W\ 

arg 1 



4T(&) 



1 



1 



4T(/3 e ) 



(77) 



Fig. 6 (a) I displays the amplitude- frequency spectrum for 
the antenna parameters given in |II A| Quantities plotted 
arc the dimensionless amplitude, A(u)/F, given by Eq. 
(77 1, versus the dimensionless frequency ui/uji^. There- 
fore, the peak of the fundamental is located at 1. The 
peak of the first collective mode (first frequency of the sec- 
ond band) is approximately located at 1^1,2/^1,1 = 119.23 
as seen on Fig. |6(b)| The dissipation parameter has been 
chosen in such a way as to provide a fundamental peak 
with quality factor Q = [it^\,\l v t on the order of ~ 10 3 , 
a typical value in experiments. The physically relevant 
part of the spectrum (corresponding to n < TV, that is 
|/3 c (w)| < (N + l/2)7r), is drawn in solid line while the ir- 
relevant part is drawn in dotted line. Because of the weak 
damping, frequency peaks occur when the driving fre- 
quency is close to one of the modal frequencies, uJ2n+i,k- 
As seen in Fig. 6(a)| peak amplitudes decrease rapidly 
as the frequency increases within a given band. But the 
first peak of the collective band, approximately located at 
u>2.i — 119.23wi.i, is much higher than the last peaks of 
the fundamental band. That is why it is easily observed 
in experiments. An analytical formula for peak heights 
is obtained in the next section in the weakly dissipative 
regime. 

Also of interest are the dips occurring slightly before 
the peaks. Their occurrence is due to the fact that at cer- 
tain frequencies, the integral of the shape of the central 
beam over its length is close to zero. This is especially 
important for detection schemes involving an electromo- 
tive voltage because the signal is close to zero in this case. 
For a non- dissipative system, there are frequencies such 
that the integral is exactly zero. The shape of the central 
beam resembles the modal shapes l2n+i(w), n > 1, in 
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FIG. 6: (a): Amplitude- Frequency spectrum 
parameters given in 



(A(uj)/F versus obtained from Eq. (77l for N 



20 and the antenna 

II A Solid lines correspond to n < N (physically relevant frequencies) and the dotted one to n > TV. (b): 
Zoom in the region of the 1st collective frequency peak. For both (a) and (b), analytical results for peaks (+) and dips (o) are 
obtained from Eqs. 1 80 1 and (811, respectively. 



this case. From Eq. (77), we see that spectral dips are (64 1 for 7 C and the result obtained in (52) yields 



given by the simple relation 4T(/3) = (3, where T{(3) is 
given in ( 70 1 . Denoting with a hat all quantities related 



to the dips, we can show that the solution of the previous 
equation reads approximately 



fhn+l — fhn+1 + 



4- A 



n> 1. 



2n+l 



(78) 



This explains why dips are close to the frequency peaks. 
The corresponding frequencies are obtained by solv- 
ing the secular equation Q(j, $2n+i) = for 7 (see 
Eq. ( pjSj l) whose solutions ^n+i.k define Cj 2n +i,k — 
\/£t/ Ht{j2n+i,k/l) 2 ■ La st, the deepest dip occurring 
around 100 in Fig. 



of the fundamental band, obtained in the non dissipative 
case for Woo,i/Vl,i = (7004/71, i) 2 — 100.103. Though, 
this part of the spectrum is not physically relevant. 



5. Weak dissipation 

Analytical results for weakly dissipative systems are 
obtained by expanding the quantities of interest around 
their value in absence of dissipation (v t = 0). Let us 
assume that, for v t = 0, ojq is the frequency of in- 
terest, a modal or dip frequency, for example. Then 
70 = lifit^o/St) 1 ^ 4 and (3 given by ( 32 1 are real. For 
uj ~ uj , we have ^(w, v t ) ~ (3q + [{to - u> ) d^/duj\ + 

v t 07c/^tU dPc/dlc\o- Now ' usin § E( l- © for E( l' 



\P»-Po\ 



St 



2Nm t 



i(7o) 



,2 , l' v tVo\ 



) 



, (uj~oj ) (79) 



Notice that, because u> ~ loq, the function under the 
square root is, in first order in v t , equivalent to (oj 2 — 
w 2 ) 2 + (vtui / /it) 2 , which is the typical form for a har- 
monic oscillator. But /3 c (w, ft) has to be expanded up to 
second order in v t and uj — uj to compute the correct fre- 
quency shift induced by the damping because the latter 



6(a) is due to the upper band edge is °f order 



By expanding (3 c (uj) around fcn+i and using 
Eq. (77), we find that the amplitude of the 
peak with frequency W2 n +i,k is given by A(u>) ~ 
16Ftanh 2 ( i 22„ + i/2)/(/3 2 „ +1 |/3 c i ( W ) - /3 2 4 „ +1 |). Therefore, 



from Eq. ( 79 ) we obtain the peak height as 



a peak 
"mi 



Vi 



1- 



2Nm t 
m b 



m.k) 



- -1 tanh 2 



(80) 

where m = In + 1, n > 0. In the same way, 
expanding /3 c (uj) around fon+i (see Eq. (78)) we 



can show after some algebra that the amplitude of 
a dip with frequency t2>2n+i,k is given by A(u>) ~ 



8 2 



i/2 



tanh(/3 2 „ +1 /2) 



1 



4 

'2n-\ 



Hence, the 
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minimum of the dip 



A 



dip 



FL A v t ( /3 m /z 



4£t Vtanh(/3 m /2) 



-1 x 



, 2Nm t r , 

1 + -L{ lm , k ) 



II 



(81) 



with m = 2n + 1 , n > 1 . 



6. Electromotive potential spectrum 



From Eqs. (75 1 and (76 1, it is now straightforward to 



evaluate the electromotive potential V em f (t) given in Eq. 
d59l. It reads 



V emf (t) _ dy(t) 



(v d A(uJd) sin(wd* - 6{u)d) - n), (82) 



from which we eventually deduce the time r.m.s. 



BFLoJd 1 
V2 \ftt d \ 



4T(/3 Cid ) 



ft, 



(83) 

From the formulas of the previous section we can evalu- 
ate the peak maxima and dip minima of the electromotive 
r.m.s. and we are able to compare them with those ob- 
served in experiments. 



IV. NONLINEAR SYSTEM DRIVEN BY TWO 
FREQUENCIES 

It has been observed experimentally [28 that, if in ad- 
dition to being driven at the frequency of the fundamental 
mode the antenna is also driven at the frequency of the 
collective mode, the frequency peak of the fundamental 
mode experiences a slight shift. This frequency shift is 
the signature of a modal coupling that occurs because of 
the presence of nonlinearity and dissipation in the sys- 
tem [27] . To explain the interaction of these two widely 
spaced modes, we investigate the effect of a two- frequency 
driving on the response of a weakly nonlinear and dissi- 
pative antenna. We supplement the equations of motion 
of our continuum model, Eqs. (18) and (19 1, with non- 



linear terms that take into account the possible material 
and geometrical nonlinearities of the structure and with 
damping terms proportional to the transverse velocity of 
the elements. These terms, small compared to the am- 
plitude of vibration of the antenna, are treated within 
the Lagrangian approach described in [27] to derive the 
frequency-amplitude relation of the model. 



A. Lagrangian approach 



(y n (x, t), i] ny k(x, £ , t)) T , to be related to two parameters: 
P n , root of cos ft cosh (3 = 1 and ^ n ,k, solution to Eq. 
(32 1 or (|35| . The frequency uj n ,k of tp n ,k is determined 
by Eq. p6l and its modal shape derived from Eqs. (29 1 



and (33). In particular, the fundamental mode is given 



by the parameters (/3i,7i,i) and the collective mode we 
are interested in by the parameters (fti,"/i l2 ). For clar- 
ity in the notations, we rename these two modes ipi — 
(yi(x,t),r)i(x,£,t)) T and <p 2 = (y 2 (x, t), r) 2 (x, £,t)) T , re- 
spectively and denote their respective frequencies by w\ 
and lj 2 where 



U), 



[St /7M\ 2 



ie {1,2}. 



(84) 



It is observed experimentally that tpi and tp 2 are coupled 
in the sense that the amplitude-frequency curve (reso- 
nance) of the fundamental mode is modified when the 
higher mode is driven. This coupling is attributed to the 
presence of nonlinearities in the antenna. To explain this 
phenomenon, we treat the nonlinearities and the damping 
affecting the system as a perturbation of the fundamen- 
tal and collective modes, tpi and (p2- These perturbative 
terms are responsible for a modulation of the linear modes 
that we evaluate by a multiple scale method. Hereafter, 
we closely follow the Lagrangian approach of Ref. [27] be- 
cause it offers a particularly suitable framework to derive 
the modulation equations. 

When the driving amplitude (or power) is small 
enough, the solution tp = (y(x,t),r](x,£ y ,t)) T to the non- 
linear equations of motion is, to a good approximation, 
given by a superposition of ipi and tp 2 with slowly modu- 
lated amplitudes. Following the multiple scale approach, 
we write it as (p = A\tp\ + A 2 (p 2 , that is 



,(5Ii))-{^ (r2 ^ (u) ^ To Uo 

+A 2 (T 2 )Y 1 (u)e^ T « f^ v) \ +c .c. 



(85) 



In this expression, c.c. denotes the complex conjugate, e 
is a small bookkeeping parameter, Hi(v) = Hi^v) where 



H is defined in Eq. (29 1 and, according to the multiple 



scale method, we have introduced two time scales, Tq = t 
and T 2 — e 2 t. The Lagrangian of our system is given by 



C 



d.r 



Mfc ( dy\ 2 £b ( d 2 y s ' 



dt 



2 \dx 2 



2N 



dx / d£ 



2 \dt dt) 2 \ae 



Li Li 

(NLT) + F l cos(ftt) J ydx + F 2 cos(wi) J ydx, (86) 



In the previous section, we solved the linear con- 
tinuum model exactly and found its modes, (p n ,k = 



where (NLT) stands for "NonLinear Terms" . To express 
the fact that the driving frequencies, and u>, are close 
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to the linear frequencies of the fundamental and excited 
modes, we write them as 



ft = u)\ + e 2 (Ji ; uj = u>2 + s a 2 . 



(87) 



To describe the nonlinear response of the cantilevers and 
the central beam, neglecting the effects of rotatory inertia 
and shear deformations, we add the following nonlinear 
terms to the Lagrangian |27) 



£bfdy_ d 2 y 
2 \ dx dx 2 




2N 



£ t ( drj d 2 r\ 



2 v<9£ oe 




Finally, we take into account the damping effects of the 
viscous forces acting on the antenna through the virtual 
work 



SW = 



dxC'y^Sy 



2N 



dx J e?£ Cjj 



dy dr] 



o o 



6r), (89) 



where C y and C n are the viscosities of the beam and the 
cantilevers, respectively. 

To account for the fact that damping effects and driv- 
ing forces are of the same order of magnitude as the non- 
linear effects, we scale the viscosities as C y — e 2 c y and 
C v = e 2 c v and the forces as Fj = e 3 fj, j = 1,2. We 
now proceed as explained in Ref . [2"T] to derive a time- 
averaged Lagrangian from Eqs. (851, (86) and (88 1. We 



substitute (851 into the Lagrangian (86 1 and also into 



the virtual work (89), perform the spatial integrations 



and keep the slowly varying terms only - i.e. those that 
are either constant or function of T 2 only. This yields: 



]T {iM jUj (A 3 A> - A'Aj) + C^Atf- 



3 = 1 



cc)}+2C 12 \A 1 \ 2 \A 2 \ 2 (90) 



(6W) 



QjSAj + cc, with Qj — 2iujjfijAj, (91) 



where 



Mi = 



t'j 



C, 



m b + 2Nm t L. h ; Tj = -^fj^A 
-{Lc y + 2Nlc n K 3J ) 



r imb oj 2 sr 2 £b 2iVL 



L 2 U> I 2 

T imb (iu 2 + w|) 6T 2 £ b 2NT 3 



L 2 



I 2 



(Kn22 + ^ 2211 + 4Xi212) 



3 33 73 ^3333 



m t (u> 2 + w|)Ji2 
(92) 



and where 



r 2 
r 3 
r 4 

Lij 
An 



[Y[(v)fdv) du, 



vo 
1 



[Y{{u)Y{'{u)Ydu, 



Y?{u)du, 



Yi(u) du, 



hi(v)hj(v) dv, 



hi(v)(hj(v) - 1) dv, 



Vo 
1 



K, 



ijkl 



h', i (y')h!^(y) dvj dv, 
K{ v ) h ' 3 {v)hl{v)h'[(v)dv. 



(93) 



with hj(v) = Hj(v) + 1. A numerical/analytical evalua- 
tion of the above quantities for the antenna dimensions 
given in |II A| is provided in appendix |C| 



B. Modulation equations 

Applying the extended Hamilton principle (see Ref. 
[27]). we obtain the equations of motion for the modu- 
lations j4i(T 2 ) and A 2 {T 2 ) as 



dC 

dT 2 KM, 



dC 

dAi 



Qi, J£ {1,2}. 



(94) 



From Eqs. (90 1 and (91), we then derive the following 



pair of modulation equations 

2*w<(M < 4 + laAi) = -2Ai (C u \Ai\ 2 + 

+ C l2 \A \ 2 ) + T t e^ T \ (95) 

where € {1,2}, i ^ j. Looking for solutions in the 
polar form 

Ai(T 2 ) = l - ai {T 2 ) e *C^-«.(T a » i e {1)2}j (96) 
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and separating the real and imaginary components of 
Eqs. (pi yields 



u l (o l - O'ijMidi = j (C lt a 2 + Ci 2 a]) - cos(6> t ), 
u i {M i a' i + Midi) = ^sin(^). (97) 

Looking for steady state (periodic) solutions, we impose 
a'j = and 9'j = an 
amplitude relations as 

1 



a'j = and O'j = and we finally obtain the frequency 



0"! 



(T 2 



4MiWi 

1 



4M 2 w 2 
together with 

tan 9 



[C lia 2 + C 12 a 2 ] ± 
[C 22 a2 + C 12 a?] ± 



J .2, ,2 Ji,.r2' 0^) 



F 2 
l 



M^a 2 M 2 



F 2 



M 2 w 2 a 2 



4 <"> 



x 2 - 
.7 



C« a 2 + C12 a 2 - AMiUjjij, 



(100) 



In first approximation the steady state solution can be 
cast into the form 



y(x,t) 



r, 1 [ i?1 1 ( w )) << ^ ( ^ / -^'>+ 



+a 2 



cos(a;t - 6> 2 ) L (101) 



The amplitudes ai and a 2 are assumed small enough for 
the perturbation expansion to hold (notice that the book- 
keeping parameter e has been absorbed in the amplitudes 
and that Eqs. ( 98 1-( 99 ) can be used as such provided the 



detunings <jj are redefined as aj = £ 2 <Jj, the viscosities 
fj,j as fij = e 2 /ij and the forces Tj as Tj = e J 7 ,-). 

C. Discussion 



The frequency-amplitude relations ( 98 1-( 99 ) allow us 



to evaluate the frequency shift of the fundamental peak 
induced by a driving of the higher mode at the exact 
linear resonance frequency, lo 2 . This frequency shift is 
determined as the difference between the maximum of the 
resonance peak of the fundamental mode and lu\. Now, 
the amplitudcai becomes maximum if the square root in 



the r.h.s. of (98 1 vanishes, that is, 



(102) 



For a system whose higher mode is driven exactly at fre- 
quency uj — u>2 we have of course tr 2 = and then the 
amplitude of the second peak is solution to 



1 

4M 2 w 2 



C 22 a 2 2 + C 



12 2 2 



± 



l 2 /7 2 

i 2 uJ 2 a 2 



(4 



= o, 

(103) 



which is a cubic equation for a 2 . Once the solution 
a 2 (^ r i,^ r 2 ) is known, we can reinstate it in (98) and we 
obtain the frequency shift, o\ — — uji as 



CTi (^1,T 2 ) 



1 



C u 



C 12 a 2 2 {T u T 2 ) 



(104) 

which provides the frequency shift as a function of the 
forces (or driving power) , T\ and Ti , of the fundamental 
and excited modes. 



V. CONCLUSION 

Here, we have presented two analytical models that 
yield closed- form solutions describing the dynamics of the 
coupled-beam resonator, dubbed the antenna structure. 
This structure is a prototype of a class of two-element me- 
chanical structures that can be envisioned in specific ap- 
plications that involve coupled mechanical oscillators and 
hierarchical structures. The inherent modifications asso- 
ciated with the dynamics of coupled-element structures 
can be engineered to result in advantageous frequency 
and amplitude performance, which is otherwise difficult 
to obtain with simple geometries. In particular, the mea- 
surements of a similar nanomechanical fabricated device 
have been reported previously |121 I23j to demonstrate 
some of the highest mechanical resonance frequencies up 
to 3 GHz, reported to date. 

The continuum model allows for a clear comparison 
of the modal shapes and spectrum with full finite ele- 
ment analysis of the structure. The findings and result- 
ing discussion of Sec. |III G| elucidate the behavior of the 
coupled-element system in the fundamental and first col- 
lective modes. In particular, it is shown that the en- 
hanced effective amplitude of the collective mode results 
from the collective excitation of the cantilever continuum 
at high frequencies, while the supporting central beam 
effectively adds the cantilever motion by moving in its 
fundamental mode shape with zero nodes, thus providing 
maximum transduction of the cantilever displacement to 
the measured magnetomotive voltage. 

We have further investigated the driven damped model 
of the system as well as the nonlinear modal coupling be- 
tween widely spaced modes of the structure using pertur- 
bation theory techniques. The results elucidate the re- 
sponse of experimentally measured structures that show 
modal interactions even in the limit of linear driving (to 
be published elsewhere). The analysis, however, is read- 
ily applied to a general set of coupled element weakly 
damped and weakly nonlinear resonators, and it illus- 
trates the nontrivial modifications in the dynamics of 
such systems that can be carefully engineered to suit spe- 
cific technical needs. As was mentioned in the section [TJ 
applications of RF MEMS and NEMS devices are numer- 
ous in the areas of wireless communications and frequency 
manipulation. 
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APPENDIX A: GENERAL SOLUTION TO THE 
DISCRETE MODEL 



The solution to the system of equations ( 12 )- (II 31) reads 



Y{u) = C[ 3) cos(fiu) + C { 2 j) sin(/m) + C { 3 j) cosh(/wt)+ 
+ C± ' sinh(^tii) for u €(%_!, Uj), (Al) 
with j e {1, . . . ,N + 1}, and 

ffi(u) = [^1(7) cos(7w)+A 2 (7) sin(7«)+yl3(7) cosh(7t))+ 
+ A 4 ( 7 )siiih(7*;)-l]l>;), (A2) 

with i € {1, . . . , iV}, ito = and ujv+i = 1- Given the 
boundary conditions ( ]14[ ) . we can solve for A; (7) in Eq. 
(A2| and find the result given in Eq. (30). Now, let 
C U) = ( C -« ) _ ; c| 3 )) T . Using the boundary conditions 



(14 1 for Y(u) at Uj, we obtain 

= (l + 2aA 2 ( 7 )K j ( M ))C« ! j e {1, 



,N}, 
(A3) 



where the quantities involved are defined in section |II C| 
The boundary conditions Eq. (fi~4]) at u — show that 

^(D ^,(D _„j „a) ^(T>^ , ■ 



-C[ 1} and CJ 



C 2 \^ The problem is then 



reduced to the two components C± and C^ 1 ■ Rewriting 
them as the 2— vector c 1 - 1 ' 



(C[ 1 \ci 1) ) T , we have 



LcW, 



(A4) 



where L is defined in Eq. (17 1. At the other end of the 
beam, u = 1, the boundary conditions can also be cast 
into a matrix form as 

T(m) C {n+1) = 0, (A5) 

where T(/i) is again defined in Eq. ( fi"7| ). Putting Eqs. 
(A3 1, (jA4h and (A5) together, we finally see that 



M(u)c (1) = 



(A6) 



where the 2x2 matrix M(w) is given in ( 16 1. This system 
has a nonzero solution cW if and only if the determinant 
of M(u>) is zero, hence the secular equation 



det T(/i) 



N 

IK 1 



2aA 2 ( 7 )Kj( M )) 



0. (A7) 



Once a solution w n of (A7) has been found, the coeffi- 
cients (w n ), i = 1, . . . , 4, of the correspondi ng m ode, 
Y n (u), are automatically determined by Eqs. (A3 1 and 
(|A6| as 



C« +1 >K) = 



JJ(l + 2aA 2 ( 7n )K k (/i„)) 



,fc=i 



c (1) K), 

(A8) 



with j e {1,.. . ,N} and 



C^kj^cOM and c (1) K)=^(_Mn^i))' 

(A9) 

where M$j are the coefficients of the matrix M. Ulti- 
mately, the normalization factor J\f is determined from 
J^Y n (u) 2 du=l. 



APPENDIX B: SOLUTION OF THE DISCRETE 
MODEL FOR N = 1 (TWO CANTILEVERS) 

In this appendix, we find the solution of the discrete 
model for two cantilevers (N =1) located on both sides 
of the beam in its middle and we briefly compare it to the 
solution of the continuum model with the same number 



of cantilevers. Using the general result provided in (A7), 
we can cast the secular equation into the form 



-4(^)5(7,/!) = 0, 



(Bl) 



where 



and 



A(n) = sin - cosh ^ - sinh - cos - (B2) 



5(7, /*) = 2aA 2 ( 7 )(cos | cosh I - 1) 



sin — cosh — h sinn — cos — . (-B3) 
2 2 2 2 v ; 



where A 2 ( , y) is given in (30 1 and a = w/W. The fac 



torization of the secular equation has a clear physical 
meaning. It is the result of the mirror symmetry of the 
problem with respect to the middle of the beam. Because 
of it, the modes of the discrete model are cither symmet- 
ric (Y s ) or antisymmetric (Y a ) with respect to the latter. 
In what follows, we use Eqs. (A9l and (Bl I to calculate 
the analytical form of the mod; 



al shape of t 



ne symmetric 



and antisymmetric modes for the central beam and the 
cantilevers. 



1. Symmetric modes 

Symmetric modes satisfy Y s (u) = Y s (l — u). Their 
secular equation is given by (/x > 0, 7 > Q), 



5(7,/i)=0, 



(B4) 



where 5(7, fi) is given in (TB3J) and, according to Eq. ( 11 ), 



fi = "fL/l. The corresponding frequency is obtained from 



(B5) 




u> = \ — V 
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The expression for the modal shape of the central beam 
is, 



Y f 



(tt) = Af( cos(fiu) — cosh(//u) + 



/ sin | + sinh £ 



V cos £ - cosh g 



x (sin(/iu) — sinh(/iM))|, 

with a normalization factor A/" given by 
6 (cos ^ cosh # - 1) 



1 



(B6) 



AT: 



/i (cos I - cosh 

x (cos — smh — h cosh — sin — 
v 2 2 2 2 ; 



-1/2 



(B7) 



To evaluate the deflection of the cantilevers, we use Eq. 
flA2"t and find 



Hi(v) = [^1(7) cos(7«) + ^2(7) sin(7?;)+ 

+ A 3 ( 7 ) cosh( 7 u) + Ai( 7 ) sinh( 7W ) - l]F s (l/2), (B8) 



where the coefficients A, are given in Eq. ( 30 1 and where 



the deflection of the middle of the beam is determined 

' 1 — cos £ cosh £ 



from ( B6 1 as 

F s (l/2) = 2Af 



COS ' 



cosh 



(B9) 



beam, symmetric modes are affected by their motion and 
the secular equation becomes 5(7, /x) = (see Eq. (B4|) 
rather t han A + (fi) = 0. Their modal shape changes from 
( 33 I to ( B6 ) . For antisymmetric modes, however, the can- 



tilevers do not move so that, modal shape and frequency 
remain unaffected. 

The first roots of Eq. (Bll I are given by fi ~ 7.8532, 
14.1371, and (2n + 1/2)tv as n becomes large. The corre- 
sponding frequencies are obtained from Eq. (B5l. Notice 



that, in this particular case, the modal shape of an an- 
tisymmetric mode of the discrete model is exactly the 
same as the modal shape of the continuum model. Their 
frequency differs however. Even though fi and j3 (see 



Eq. (34 1) satisfy the same secular equation (and are thus 
equal), the frequency of the continuum model is obtained 
by solving equation (32 1 whose solution, 7 C , is not pro- 



portional to (3 while the solution of the discrete model, 7^, 
satisfies 7t ; = fil/L. The frequencies, both obtained from 
Eq. (B5l, are thus different. The exact similarity of the 



modes may seem surprising at first glance because, if the 
cantilevers are at rest in the discrete model and thus do 
not participate to the motion, in the continuum model, 
the force density of the cantilever has been spread all 
over the beam and, consequently, the cantilever contin- 
uum moves with the beam. Nonetheless, as the force den- 
sity it generates is everywhere proportional to the mode 
shape Y(u), the frequency is detuned but the mode shape 
remains as is. 



2. Antisymmetric modes 



3. Comparison with the continuum model 



Antisymmetric modes satisfy Y a (l — u) = —Y a (u), 
which leads to Y a (l/2) = 0. The deflection of the can- 
tilevers is then zero, H±(v) = 0. In this case, we di- 



rectly see from equations (12 1 and (14), that Y a (u) sat- 



isfies a simple clamped-clamped beam equation and that 
its third derivative has no discontinuity in u — 1/2, given 
that H\(v) = 0. The shape of the antisymmetric modes 
is then similar to Eq. ( 33 ) , which yields 



Y a (u) — I cos(/iu) — cosh(/i«) 



cos fi — cosh fi 



sin fi — sinh fi 
(sin(/iu) - sinh(/iu))|, u E [0, 1]. (B10) 



From Eq. (Bll, we see that [i has to be a root of the 



secular equation 



A(fx) = 0, 



(Bll) 



where A{jj) is given in (B2|. It turns out that Eq. (BllJ 



is equivalent to the usual clamped-clamped secular equa- 
tion, cos(/i) cosh(/i) = 1, provided it is restricted to anti- 
symmetric modes. Indeed, cos(^) cosh(/i) = 1 factorizes 
as A(fi)A + ([i) = where A + ([i) = sin fi/2 cosh \ij2 + 
sinh n/2 cos /x/2: A(fi) — provides the frequencies of 
antisymmetric modes while A + (/i) = gives the frequen- 
cies of the symmetric modes of a simple clamped-clamped 
beam. When cantilevers are affixed to the middle of the 



According to the dimensions of the antenna, L/l ~ 21.4 
and w/W ~ 0.5. Then fi ~ 21.47 and a = 0.5. Plug- 
ging these values into Eq. (B4|, we can solve for 7^ and 



obtain the values for the symmetric modes reported in 
Table [I] namely j dA = 0.2149 and j d>3 = 0.5035. For the 
antisymmetric modes, we use fi ~ 21 Aj and the values 
quoted in the previous section and find 7^2 = 0.3670 
and 744 = 0.6606. Now, solving equation (34 1 for 
f3 = f3\ to /?4, where j3 n are the successive roots of 
cos(/3) cosh(/3) = 1, we find the corresponding values of 7 C 
for the continuum model, reported in the third column of 
Table [T] As we see, these values are very close to those of 
the discrete model. Using Eq. (B5l we finally report the 



frequencies of the discrete (4th column) and continuum 
(5th column) modes. Their relative frequency difference 
is displayed in the 6th column. The modal shapes of the 
beam and the cantilevers from the fundamental mode up 
to the third excited mode are depicted in Fig. [7J For the 
cantilever continuum, we have displayed the deflection at 
u = 1/2 only. As we can see, symmetric modes are very 
similar and antisymmetric modes are exactly the same. 

Beside the first levels presented in table [i] we have 
found numerically that the agreement between the spec- 
tra of the discrete and continuum models is very good 
in the sense that there is almost always one frequency of 
the continuum spectrum that closely matches the corre- 
sponding frequency of the discrete model. The shape of 
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TABLE I: Values for the parameters jd (discrete model) and 7 C (continuum model) derived from Eqs. ( B4 1-( Bll I and { 34 1 , 
respectively. Frequencies are determined from Eq. ( B5 1 for the antenna parameters given in the introduction. The symmetry 
class of the modes is denoted by s (symmetric) or a (antisymmetric). 



Mode 


Symmetry 


Id 


7c 


fd a (MHz) 


fc a (MHz) 


A/// (in %) 


Fundamental 


s 


0.2149 


0.2185 


31.81 


32.88 


-3.25 


First Excited 


a 


0.3670 


0.3628 


92.74 


90.64 


2.31 


Second Excited 


s 


0.5035 


0.5079 


174.6 


177.7 


-1.74 


Third Excited 


a 


0.6606 


0.6530 


300.5 


293.4 


2.33 



"The frequency is determined as / = u)/2ir 




(a) (b) (c) (d) 




FIG. 7: Fundamental to third excited mode shapes of the central beam ((a),(b),(c),(d)) and of the cantilever ((e),(f),(g),(h)). 
The results of the discrete and continuum models are in red and green, respectively. For comparison, the deflection of the 
cantilever continuum is shown at u = 1/2 only. The antisymmetric discrete and continuum modes are exactly the same but 
their frequencies are different (see Table IB. 



the central beam is very similar (same n) for matching 
frequencies. The frequencies of the discrete model, how- 
ever, are quite regularly spaced while in addition to pro- 
ducing these regularly spaced frequencies, the spectrum 
of the continuum model also forms clusters of nearly de- 
generate frequencies as explained in section |IIIE| It is 
clear from the data we have obtained that, in the dis- 
crete model, the excitation level (n) of the central beam 
increases regularly with the frequency. In the continuum 
model, however, this excitation level increases from n — 1 
to infinity within each band, that is from toi.k up to Woo,fc 
determined by Eq. (38 1. Therefore, the continuum fre- 
quencies that match the discrete ones are located away 
from the band edge clusters of the continuum spectrum. 
The continuum frequencies that accumulate close to the 
band edges have therefore no (N = 1) discrete equivalent. 
In other words, collective modes such as those observed 



in the continuum model appear only when the number of 
cantilevers is large enough for the continuum approxima- 
tion to hold. 



APPENDIX C: ANALYTICAL AND NUMERICAL 
RESULTS FOR EQS. ([93]) 



Using Eqs. (33) and (29 1, all the integrals involved 



in the calculation of the effective parameters of the time 
average Lagrangian, see Eq. d93j), can normally be evalu- 
ated analytically. For I\, i £ {Tf - • • , 4}, results are simple 
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enough to be displayed below 
Pit 



r 2 
r 3 
r 4 



2 
fit 
10 
3 
4 
4f 

Pi 



(Pit + 2) ~ 6.1513, 
(5/?i* + ll) ~ 2846.4975, 

Of 3 

3-t 4 - — ) - l.s.VH). 



= — ~ -0.8308 



where t — tan(/3i/2). Note that cos(/?i) = l/cosh(/3i) 
and sin(/3i) = — tanh(/?i). Apart from = La — 
2A 2 (7 l )/7 l and L« = £(7i,i) ; given in Eq. ( |5"lj l, the 
other integrals have been evaluated numerically for the 



parameters given in II A We have found 



L n ~ 1.00012 ; L 22 ~ 3.89887, 

A n ~ 6.16 x 10" 5 ; A 22 ~ 4.9687, 

J n ~ 1.64 x 10~ 16 ; I 22 ~ 232.49, 

I 12 ~ 1.94 x 10~ 7 ; X1212 ^ 8.26 x 10~ 7 , 



i^iiii 

^1122 



6.58 x 10~ 16 ; ^2222 
9.84 x 10~ 7 ; K 22n 2 



± 1.07 x 10'\ 
6.98 x 10~ 7 . 



So that, finally, 



Mi = 1.74 x 10~ 14 ; 
Ti = -4.44 x lO" 6 ^ 
C*n = -5.07 x 10 13 
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M 2 = 4.17 x 10" 14 , 
; C12 = 1.65 x 10 17 , 
C 22 = 1-05 x 10 21 , 



wi = 1.55 x 10° ; cu 2 = 1.85 x 10 



10 



where all the results are given in SI units. From the quan- 
tities u>j, the frequencies of the fundamental and collec- 
tive modes are 24.7 MHz and 2.94 GHz, respectively. 
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Because of Eqs. (|A2} and (30), we have Y"'(uj) 
Y'"(uJ) = -RH?'(0) = 2 1 3 RA 2 ( 7 )Y(u J ) 
A/j, 3 aA2(-f)Y(uj), where a — w/W. 



